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A GENEALOGY OF CONVEX SOLIDS VIA LOCAL AND 
GLOBAL BIFURCATIONS OF GRADIENT VECTOR FIELDS 


GABOR DOMOKOS, PHILIP HOLMES AND ZSOLT LANGI 


Abstract. Three-dimensional convex bodies can be classified in terms of the 
number and stability types of critical points on which they can balance at rest 
on a horizontal plane. For typical bodies these are nondegenerate maxima, 
minima, and saddle-points, the numbers of which provide a primary classi¬ 
fication. Secondary and tertiary classifications use graphs to describe orbits 
connecting these critical points in the gradient vector field associated with each 
body. In previous work it was shown that these classifications are complete in 
that no class is empty. Here we construct 1- and 2-parameter families of con¬ 
vex bodies connecting members of adjacent primary and secondary classes and 
show that transitions between them can be realized by codimension 1 saddle- 
node and saddle-saddle (heteroclinic) bifurcations in the gradient vector fields. 
Our results indicate that all combinatorially possible transitions can be real¬ 
ized in physical shape evolution processes, e.g. by abrasion of sedimentary 
particles. 
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1. Introduction 

1.1. Motivation and background. The evolution of shapes of abrading bodies, 
such as pebbles in river beds and on beaches, has been studied for over 70 years 
(e.g. 23; [24l t25j 12| 13]). Data from NASA’s Curiosity Rover on Mars [29, [3j| 

has rekindled interest in the subject. In addition to classical shape indices such as 
axis ratios and roundness mm, a recent approach considers the evolution of the 
number of static equilibrium points N(t ) on the surface of an abrading body, i.e., 
points on which the body can balance at rest on a horizontal plane [Ml [9J1 [5]. 
Unlike shape indices, which require length measurements, the integer N(t) can be 
counted in simple experiments [9]. 

Abrasion occurs primarily on a body’s convex hull, so to formulate a precise 
and relatively simple model we restrict our analysis to convex bodies K of uniform 
density, with surfaces described by scalar Euclidean distance functions vk measured 
from the center of mass Ck- For such bodies static equilibria are critical points of 
tk at which the gradient Vr^ = 0. 

The surface dK of a generic convex body K can exhibit three types of nonde¬ 
generate critical points: local minima, maxima and saddle-points, which are sinks, 
sources and saddles of the gradient vector field v = Vr^- Let S', [/, H respectively 
denote the number of each of these points. Since dK is a topological 2-sphere, the 
Poincare-Hopf Theorem [2] implies that 

(1) S + U-H = 2. 




FIGURE 1. Primary equilibrium classes. Left: examples of convex 
bodies; rows and columns correspond to the numbers S and U of 
sinks and sources, respectively. Right: the ‘Columbus Algorithm’ 
of |28j defines a hierarchy among primary classes. Arrows indi¬ 
cate arbitrarily small truncations of the convex body, creating one 
additional sink or source and a saddle-point. 

The classification schemes introduced in [28! are based on these numbers. Specif¬ 
ically, the primary class of a generic convex body K is defined as the pair of integers 
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{5, U}. In [28; it was shown that no primary class {i,j} is empty and a hierarchy 
among these classes was defined via the Columbus algorithm. Using explicit trun¬ 
cations that remove small portions from K by slicing along convex surfaces, this 
algorithm generates a pair of convex bodies K' G {i + 1, j} and K" G {z, j + 1}, 
as shown in Figure [l] Thus, starting from the gomboc {1,1}, every row and col¬ 
umn can be populated, implying that the primary classification is complete in this 
‘static’ sense. 



Figure 2. Graph representations of the gradient flow on the tri- 
axial ellipsoid in primary equilibrium class {2,2}. (a) Distance 

function tk given in spherical polar coordinates, (b) 3-colored 
quadrangulated primary representation Q 3 (v). (c) 3-colored tri¬ 
angulated representation T 3 (v). (d) Quasi-dual, 2-colored quad¬ 
rangulated representation Q 2 (v). The colors refer to vertices, iden¬ 
tifying them as sinks, sources, and (in (b,c)) saddles. 

More refined methods exist for classifying the properties of gradient vector fields 
v = Vr^, including graph representations of their Morse-Smale complexes m- 
The vertices of these graphs are fixed points of v, and the edges can be either iso¬ 
lated heteroclinic orbits connecting saddle-points, representative non-isolated het¬ 
eroclinic orbits connecting saddles with sinks and sources, or both. These are 
called, respectively, the primary representation Q 3 (v), the triangulated represen¬ 
tation T 3 (v), and the quasi-dual representation Q 2 (v) Figure [2] illustrates these 
representations for the tri-axial ellipsoid. For brevity, we call all three types the 
topology graphs associated with v. Note that all three graphs Q 3 (v),T 3 (v), and 
Q 2 (v) are embedded on § 2 and we will also consider their abstract, non-embedded 
versions Q 3 (v),T 3 (v) and Q 2 (v). We remark that an abstract graph may have 
several, non-homeomorphic embeddings in § 2 . Precise definitions will be given in 
Section [2j and these graphs will play a key role in the paper. 

We call the class of convex bodies with isomorphic abstract graphs the secondary 
equilibrium class and the class of convex bodies with homeomorphic embedded 
graphs the tertiary equilibrium class associated with K. See Figure [3][a), which 
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also illustrates that a primary class can contain different secondary classes: e.g., 
the ellipsoid is not alone in class {2, 2}. In [8] it was shown that the secondary and 
tertiary schemes are also complete in the sense that no secondary or tertiary class 
is empty. 
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Figure 3. Secondary and tertiary equilibrium classes, (a) Sec¬ 
ondary and tertiary classes are contained in primary classes, (b) 
Metagraph Q with vertices at tertiary classes and edges correspond¬ 
ing to codimension 1 bifurcations; thin edges: saddle-nodes, thick 
edges: saddle-saddle bifurcations. Note that all illustrated sec¬ 
ondary classes contain one tertiary class, so in the figure the ver¬ 
tices of the metagraph Q correspond simultaneously to secondary 
and to tertiary classes, (c) Quasi-dual topology graphs Q 2 (v) of 
the tertiary classes labeled a through k in panels (a,b) and illus¬ 
trations of a through g as convex bodies 


One can ask whether transitions between different primary, secondary and ter¬ 
tiary classes are possible within generic families K( A) of smooth convex bodies, 
parametrized by A, as their shapes change. In generic one-parameter families of 
gradient vector fields only two codimension 1 bifurcations occur: saddle nodes and 
saddle-saddle connections, and they do so at isolated, critical values A = A? r M- 
Saddle-nodes involve local changes in topology in which pairs of non-degenerate 
equilibria, either a saddle and a sink or a saddle and a source, emerge or disappear. 
Saddle-saddle bifurcations are global bifurcations at which an orbit connecting two 
saddle-points exists, but the numbers and types of equilibria do not change. In the 
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former, one of the integers S', U characterizing the primary class of K increases or 
decreases by one; in the latter, the primary class remains unchanged. 



(a) 




Figure 4. Topology graphs (a) and (c) corresponding to two ver¬ 
tices of Q in primary equilibrium class {4,4} and the tertiary edge 
(b) connecting them. Graphs are shown in the triangulated repre¬ 
sentation of graph class T 3 . Note the saddle-saddle connection on 
(b), and that (a) and (c) are isomorphic as abstract graphs, but 
not homemorphic as embedded graphs on § 2 . 


To visualize these transitions we introduce the metagraph Q with vertices rep¬ 
resenting the embedded topology graphs associated with generic gradient fields on 
§ 2 . The metagraph distinguishes primary edges , on which saddle-node bifurcations 
occur between vertices in different primary classes, from secondary edges that con¬ 
tain saddle-saddle bifurcations between vertices within the same primary class but 
in distinct secondary classes, and tertiary edges between vertices in the same sec¬ 
ondary class. In Figure |3jb) primary and secondary edges are identified by thin and 
thick lines respectively. Figure [3] shows only primary classes with values S + U <6; 
here tertiary edges cannot be illustrated since secondary classes with so few critical 
points contain only one tertiary class, and therefore have unique embeddings on the 
2-sphere. Figure [4] illustrates a tertiary edge conecting two vertices in the primary 
class {4,4}. 

Our goal, which we will formally define in subsection |1.2[ is to show that not 
only the vertices, but also the primary and secondary edges of Q can be represented 
by convex bodies, i.e., physical processes describing their shape evolution may be 
represented by paths on Q. It is an intriguing question which of these paths are 
preferred by physical abrasion processes. The investigation of this question lies 
beyond the scope of the current paper, but we have some physical intuition on how 
primary classifications of convex bodies might evolve. In |[5j it is shown that the 
evolution of S and V under the partial differential equations governing collisional 
abrasion processes can be modeled by letting S and U be random variables whose 
expected values decrease with time. While this trend has been verified both in 
laboratory experiments [6] and in the field m, almost no pebbles in the primary 
classes 1}, (i,j = 1,2,...) have been found. In [7j a purely geometrical 

reason for this phenomenon was pointed out. The difficulty of reducing either the 
number of sinks or sources can be measured by the fraction of volume that must 
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be removed from a solid, referred to as robustness. In [7] it was shown that the 
robustness of classes {2, i}, {j, 2} is maximal so it is very unlikely that any natural 
pebble in these classes will be transformed into any of the classes {1 , i}, {j, 1} by 
natural abrasion. 

Much less is known about the secondary, let alone tertiary classification of con¬ 
vex bodies, and theories for their evolution by abrasion are lacking. Nevertheless, 
as for the primary case, field data indicate that natural shapes strongly prefer some 
secondary classes while other classes remain virtually empty. Already in m it was 
observed that coastal pebbles tend to be ellipsoidal. While Rayleigh [23, 24] [25j 
Bloore [3] and Firey [T2 ultimately showed that the classical exact ellipsoid is 
not an attracting state in collisional abrasion processes, nearly ellipsoidal shapes 
nonetheless dominate pebble beaches. Without exception, all those shapes in pri¬ 
mary class {2,2} for which the secondary classes were determined have topology 
graph ‘d’ of Figure [3^c) , while the other secondary class ‘c’ in {2, 2} appears to be 
missing. Similar observations apply to other primary classes. 

As a first step towards understanding these phenomena we show that the sec¬ 
ondary classification scheme of [8J is also complete in the following ‘dynamical’ 
sense. Primary and secondary edges in the metagraph £/, containing codimension 
1 saddle node and saddle-saddle bifurcations respectively, exist in the space of gra¬ 
dient vector fields v = Vvk on the 2-sphere associated with convex bodies K. 
In the next subsection we use the metagraph Q to formulate our statements and 
relate them to earlier results. Before doing so, we note that the gradient vector 
field v = \7vk cannot describe the Newtonian dynamics of the body K rocking 
on a horizontal plane, which would require a system of second order differential 
equations, but that the stability types of its fixed points correctly reflect those of 
the equilibria of K. 

1.2. Definitions and main result. We first define the metagraph Q , whose ver¬ 
tices are embedded topology graphs representing tertiary equilibrium classes as¬ 
sociated with the Morse-Smale complexes [20 of gradient vector fields of convex 
bodies. For simplicity we use the primary representation Q 3 (v), but the triangu¬ 
lated or quasi-dual representations may also be used to construct Q. The edges of Q 
correspond to codimension 1 bifurcations connecting these classes, and all possible 
one-parameter families of gradient vector fields on the 2-sphere appear in Q. We 
define the edges and vertices of Q and we will use these concepts to formulate our 
results and relate them to earlier results. 

Definition 1. Two vector fields v and w on the 2-sphere are topologically equivalent 
MM if their embedded topology graphs (of the same type) are homeomorphic. 

As noted earlier, in case of generic vector fields [26], topology graphs can be 
defined by the Morse-Smale complex associated with the vector field m- Now we 
proceed to define the metagraph Q. 

Definition 2. A vertex of Q is an embedded topology graph Q 3 (v) on § 2 , associated 
with the Morse-Smale complex of a generic gradient vector field v on § 2 . 

Definition 3. The primary class of a vertex Q 3 (v) is the pair of integers {S', I/}, 
where S and U denote the number of sinks and saddles of v. The secondary and 
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tertiary class of a vertex are the abstract graph Q 3 (v) and the embedded graph 
Q 3 (v), respectively, both associated with the Morse-Smale complex of v. 

Definition 4. An edge of Q is a one-parameter family v(A), A E [0, 1] of gradient 
vector fields connecting two distinct vertices Q 3 (v(0)) and Q 3 (v(l)) of Q. We 
require that v is generic except for a unique value A = A* E (0,1), for which v(A*) 
exhibits a codimension 1 bifurcation M- 

Definition 5. We call an edge v(A),A E [0,1] primary if the primary classes of 
Q 3 (v(0)) and Q 3 (v(l)) are different. We call an edge v(A),A E [0,1] secondary 
if the primary class of Q 3 (v(0)) and Q 3 (v(l)) are identical, but their secondary 
classes are different. We call an edge v(A), A E [0,1] tertiary if both the primary 
and the secondary class of Q 3 (v(0)) and Q 3 (v(l)) are identical. 

Definition 6. We call a vertex Q 3 (v) of Q physical if there exists a convex body 
K such that Vr^ is topologically equivalent to v. 

Definition 7. We call a primary, secondary or tertiary equilibrium class physical 
if it contains at least one physical vertex. 

Definition 8. We call an edge v(A), A E [0, 1] of Q physical if there exists a one- 
parameter family K( A), A E [0,1] of convex bodies such that Vr^(A) is topologically 
equivalent to v(A) for all values of A E [0,1]. 

Now we can formulate earlier and current results. Regarding primary equilibrium 
classes we have 

Theorem 1. All primary classes are physical 

This result, proved in [28], was generalized in [8] to include secondary and tertiary 
classes: 

Theorem 2. All vertices of Q are physical. 

In the current paper our goal is to further extend Theorems [l] and [2] by proving 
the physical existence of an important subset of edges of Q\ 

Theorem 3. All primary and secondary edges of Q are physical 

1.3. Sketch of proof. As noted above, the local truncations constructed in [8] 
modify the Morse-Smale complex of K to produce one-parameter families of con¬ 
vex bodies in which either S or V is increased by 1. However, these families do 
not (necessarily) represent edges in the metagraph Q since the genericity of the 
bifurcation was not guaranteed by the construction in [81. One-parameter families 
connecting vertices at the ends of secondary edges (saddle-saddle bifurcations) were 
not even discussed in [8]. 

Here we extend these results by constructing a 2-parameter family of convex bod¬ 
ies whose gradient vector fields are generic in the sense that certain codimension 1 
subsets (curves) in the parameter plane correspond to vector fields with codimen¬ 
sion 1 local saddle-node and global saddle-saddle bifurcations, forming primary and 
secondary edges of Q. We also show that the codimension 1 bifurcation curves meet 
in a codimension 2 saddle to saddle-node bifurcation point. 

Because secondary edges correspond to codimension 1 global saddle-saddle bi¬ 
furcations, the local methods of [8] do not apply directly. Rather, we achieve our 
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goal in two steps. In Section [ 2 ] we prove Combinatorial Lemma [lj stating that any 
secondary edge of the metagraph Q bounds a triangular face of Q of which the two 
other edges are primary. As shown in Figure [9] below, the vertices of such a face 
represent three topology graphs that lie in adjacent primary classes of Q. The tri¬ 
angles (b,c,d) and (f,j,k) in Figure[3|b) above provide examples. We then appeal to 
dynamical systems theory m in Section[3]to show that such a triangular face could 
contain a codimension 2 bifurcation point for the gradient flow v = Vr^( a) and 
describe how codimension 1 saddle-node and saddle-saddle bifurcations emanate 
from this point. 


In Section [4] we take the second step, providing an explicit geometrical construc¬ 
tion that realizes the codimension 2 bifurcation via an arbitrarily small truncation 
of K depending on two parameters. First, in Subsection |4.1| we prove that a trun¬ 
cation exists under the assumption that the resulting displacement of the body’s 
mass center has no effect on the topology of its gradient flow. Then, in Subsection 


4.2 we construct a simultaneous , auxiliary truncation such that the mass center 


remains fixed under the combined truncations, implying that the topology of the 
flow is preserved. Finally, in Section [5] we summarize our results and point out 
some possible consequences. 


2. Combinatorial part 

Before stating the combinatorial lemma, we define three classes of graphs associ¬ 
ated to Morse-Smale complexes on the 2 -sphere, of which the graph representations 
introduced above and illustrated in Figure [ 2 |b-d) are examples. As in [ 8 ], we denote 
by a quadrangulation a finite planar undirected multigraph on the 2 -sphere in which 
each face is bounded by a closed walk of length 4 (cf. [US])- A multigraph contains 
no loops but may have multiple (parallel) edges, and it is usually permitted that the 
boundary of a face may contain a vertex or an edge of the graph more than once 
(e.g. the faces with saddle-source and source-sink connections in Figure HM). 
In addition, we follow Archdeacon et al. pQ and regard the path graphs (cf. [13]) 
P 2 and P 3 as quadrangulations, where P& denotes a tree on k vertices, each with 
degree at most 2 . 

Dong et al. m introduced three different kinds of graph to represent a Morse- 
Smale complex on the 2 -sphere, as follows: 

• Q 2 is the class of 2 -vertex-colored quadr angulations. Note that as no quad¬ 
rangulation contains odd cycles, each is 2 -colorable (cf. [TJ[2T]). Further¬ 
more, the coloring of the graph is unique up to switching the colors. 

• Q 3 is the class of 3-vertex-colored quadr angulations with deg(p) = 4 for 
any p G P, and \S\ + \U\ — \H\ = 2 , where <S, U and H denote the sets of 
vertices of each given color. 

• T 3 is the class of 3-vertex-colored triangulations with deg(p) = 4 for any 
p G P, and \S\ + \U\ — \H\ = 2 , where <S, U and H denote the sets of vertices 
of each given color. 

Examples of each class appear in Figure [5j panels (c,a,b) respectively. 

It was shown in [IT] (cf. [30]) that a Morse-Smale complex on the 2 -sphere can 
be uniquely represented by a 3-vertex-colored quadrangulation in Q 3 , where the 
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• sink 
o source 
© saddle 


a) 


b) 


c) 


Figure 5. Different representations of a part of a Morse-Smale 
complex, (a) Primary topology graph in class Q 3 . (b) Triangulated 
topology graph in class T 3 . (c) Quasi-dual topology graph in class 

Q 2 . 


vertex colors represent the 3 types of critical points (maxima, minima and saddles) 
and edges correspond to stable and unstable manifolds: isolated integral curves 
that end and start at saddle points. Each quadrangle is bounded by a closed walk 
consisting of a source, a saddle, a sink and a saddle in cyclic order around the face 
and, and every saddle has degree 4; see Figure |6ja). Following Dong et al. [lOj we 
call this the primary topology graph. 

Saddle points can be removed from the primary graph without losing informa¬ 
tion [10]: first we connect sources and sinks inside each quadrangle, producing a 
triangulated topology graph in class T 3 ; we then remove all saddle points and edges 
incident to them, as in Figures |5jb,c). Since non-degenerate saddles have degree 
4, the resulting graph is a 2-vertex-colored quadrangulation in class Q 2 : the quasi¬ 
dual topology graph (cf. pH]). Here we use the latter; however, in Section [4j the 
primary graph representation is preferable. All three representations are equivalent 
in the sense that they are mutually uniquely identified. 

Let F = (pi,P2 ?P3 jP4 ) be a face of any Q G Q 2 (cf. Figure [6|a), left). Pairs 
of vertices, and/or edges connecting them, may coincide. Nonetheless, a quasi¬ 
dual representation admits only two kinds of coincidences: two diagonally opposite 
vertices, say P 2 and p 4 may coincide, and in this case two consecutive edges, say 
(p4,Pi) and (pi,P2) may coincide: these two cases are illustrated in Figures [6]^b) 
and (c), left. Note that in Figure |6jc) the internal domain bordered by the edges 
(p4,Pi) and (pi,P2) is not a quadrangular face and necessarily contains at least one 
additional vertex, as indicated by the triangles in the inner region. Figure [6^d) 
shows the remaining two exceptional cases: the trees Q = P 3 and Q = P 2 . 

The algorithm in [8] is based on repeated application of a combinatorial graph 
operation called face contraction (cf. [TJ [4] or [22]). Applied to the face F de¬ 
fined in the previous paragraph, this operation results in the contraction of the 
vertices p\ and ps into the same vertex, and the disappearance of F\ the mod¬ 
ified graphs, depending on the ‘shape’ of the original face F are shown on the 
right of each panel in Figure [6] The inverse operation of face contractions is called 
vertex splitting. Combinatorially, for graphs with at least three vertices it can 
be defined as follows. Let p be a vertex of the quadrangulation Q, with adja¬ 
cent edges Ei, £ 2 , • • •, FT, T/c+i = -ET in counterclockwise order, and note that the 
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other endpoints of some of these edges may coincide. Choose two, not necessar¬ 
ily distinct edges: E x and E y . Then we split p into two vertices pi and £> 3 , and 
E x and E y into two pairs of edges E x ,\ and E x ^ 3 , and E y% \ and E y , 3 , such that 
E Xi i, E x +i, • • •, Ey-^Ey^ are connected to p u and E V:3 , E y+1 ,..., E x _ 1 and E Xi 3 
are connected to p%. This operation can be naturally modified for primary and tri¬ 
angulated representations: in the primary representation, instead of two edges we 
choose two (not necessarily distinct) faces, whereas in a triangulated representation 
we choose two edges connecting a sink and a source. 



lPl=P3 

m=P4 

• Pl=P3 
6 Pi = Pa 


Figure 6. Face contractions on a graph in class Q 2 . As in [4], 
triangles incident to some vertices indicate that one or more edges 
may occur at that position around the vertex. Here p\ and p% are 
sinks; analogous face contractions, each removing a source from 
the graph, can be performed by switching the colors. 


In a quasi-dual graph, a transition via a saddle-saddle connection can be realized 
as a diagonal slide [ 22 ], defined as follows: consider two faces (tq, ^2,^3, ^4) and 
( 114 , 113 , 115 , H 6 ) of the quadrangulation sharing an edge (113, tq), and replace this edge 
either by (111,115) or (v 2 ,vq). Then the two faces (si, 1/3, S2, 1 ^ 2 ) and (111,53,1/2,52) 
are replaced by (si, 1x3, 53 , 1x2) and (1x1, 53 , 1x3, 52 ). 

To formulate the lemma, we need the following definition. 

Definition 9. Let Q G Q 2 be a quadrangulation. Two vertex splittings W and 
W' of Q are called twin , if: 

• The same vertex p is split. 

• Let Ai and A 2 denote the sets of edges connected to the two split vertices 
in IF, and define A[ and A' 2 similarly for IF'. Then A\ differs in exactly 
one element from A[ or A' 2 . 

In this case Q is called the ancestor of the two split graphs. This definition can 
also be naturally interpreted for primary and triangulated representations. 

Note that the second property in Definition [9] implies that A 2 also differs in 
exactly one element from A[ or A! 2 . The graphs in columns B and C of rows (a) 
and (b) in Figure [7] can be obtained from the graph in column A of the same row 
via twin vertex splittings, but the graphs in columns B and C of Figure HD are 
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isomorphic and hence have no ancestor graph. This corresponds to the degenerate 
case 3 in the Proof of the Combinatorial Lemma Q] below. 

Lemma 1 (Combinatorial Lemma). Let B,C G Q 2 be embeddings of two distinct 
abstract graphs B,C in § 2 ; respectively, such that there is a diagonal slide that 
transforms B into C . Then there is an embedding A E Q 2 and a pair of twin 
vertex splittings Wb and Wc of A such that Wb transforms A into B, and Wc 
transforms A into C . 

We remark that, as we will see in the proof of Lemma [l] there are diagonal slides 
between non-homeomorphic drawings of the same graph which cannot be derived 
from the same ancestor via twin vertex splittings. Note also the essential condition 
that the abstract graphs 5, C should be non-isomorphic; this condition excludes 
tertiary edges from our argument. 


Proof To simplify the proof, we use the triangulated variants of B and C, which 
with a little abuse of notation, we also denote by B and C. Let the two saddles 
that are connected by the saddle-saddle bifurcation be denoted by hi and ft^- This 
edge belongs to two faces of L>, say (s2,ixi,/ii) and (s2, ixi, ft-2), and similarly, two 
faces of (7, say (si,U 2 ,hi) and (si, 1 / 2 ? ^ 2 )- We note that, due to the degeneracy 
of the graph, some of the vertices or edges may coincide; nevertheless, due to the 
saddle-saddle connection, hi and ft ,2 are distinct. 



Figure 7. The connection between diagonal slides and twin ver¬ 
tex splittings for T 3 -class graphs. The graphs in row (c), columns 
B and C are isomorphic, and thus have no common ancestor, hence 
no graph appears in column A of row (c). See the proof of Lemma[l] 
for further details. 

Case i, si and S 2 , and also u\ and u 2 are distinct. Figure [7] row (a) shows 
the corresponding faces of B , the saddle-saddle bifurcation, (7, and the common 
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ancestor A from left to right. Face contractions are carried out by collapsing the 
edges (si, h i) and (s 2, h 1) into a single vertex 8, and the dotted edges starting at 81, 
hi and 82 are contracted into the single dotted edge of A. Furthermore, the edges 
(<§2, xxi), {hi,ui) and (si, u\) are contracted to (s, ui) in F?, whereas (82,1/2), (/ii, ^2) 
and (81,1x2) are contracted to (8,1x2) in C. Since (8,1x1) and (8,1x2) are consecutive 
edges of A in the quasi-dual representation, the vertex splittings belonging to the 
two face contractions are indeed twin. We remark that in Case 1 another ancestor 
can be found by contracting (1x1, /12) and (1x2, ft^)- 

Case 2 , exactly one of the pairs {81,82} or {1x1,1x2} coincide. Without loss of 
generality, we may assume that 1x1 = 1x2 = ix and 81 7^ 82. Note that as the degree 
of a saddle point is 4 , in this case there are two edges starting at h 2 and ending at ix. 
Figure [ 7 | row (b) shows the corresponding faces of F?, the saddle-saddle bifurcation, 
C, and the common ancestor A from left to right. Face contraction is carried out 
by collapsing the edges (81, hi) and (82, hi) into a single vertex 8. 

Case 3 , 81 = 82 = 8, and 1x1 = 1x2 = ix. In this case B and C are isomorphic 
graphs: Figure [ 7 ] row (c). We note that in this case the two edges starting at 8 and 
ending at ix may also coincide. 

□ 


3. Dynamical part 


In this section we describe how codimension 1 saddle-node and saddle-saddle 
bifurcations can meet in a codimension 2 bifurcation of a gradient vector field v 
on § 2 . Such a bifurcation point can be associated with each triangular face of the 
metagraph Q having two primary edges and one secondary edge. We construct an 
explicit polynomial function (x, 1/), depending on two parameters /Xi,/X2, that 

captures the behavior of v near a degenerate saddle-node whose strong stable man¬ 
ifold contains one branch of the unstable manifold of a non-degenerate (hyperbolic) 
saddle point. The parameters /xi, /X2 provide local coordinates on the face of the 
metagraph near the codimension 2 point. Since the saddle-saddle or heteroclinic 
connection is a global phenomenon, our vector field will necessarily be non-local, but 
we can nonetheless find a cubic potential function that captures the local saddle- 
node and the global heteroclinic connection. A homoclinic orbit to a saddle-node 
bifurcation point was previously shown to occur in the averaged equations for the 
periodically forced van der Pol oscillator [ 16 ], cf. [ 14 ] Sec. 2.1,Figs. 2 . 1 . 2 - 3 ]. 


We first recall the normal form of an isolated codimension 1 saddle-node in a 
gradient vector field on the plane, which can be described by a potential function 
depending on one parameter 0: 


( 2 ) V^(x,y) = y + y -Vix, 

The corresponding vector field 


( 3 ) 



gVjfl 

dy 


= -X 2 + /il, 

= -y, 


has no fixed points for /xi < 0, a saddle-node at (x, y) = (0, 0) for /xi = 0, and a 
hyperbolic saddle and a sink at (x,y) = (—^//xT, 0) and (+-071,0) respectively for 

/il > 0 . 




GENEALOGY OF CONVEX SOLIDS 


13 


We now add further cubic terms and a linear term containing another parameter 
li 2 to V^ to produce a second hyperbolic saddle that can be displaced relative to 
the saddle and sink described above. We set 

3 9 3 

( 4 ) V H1,H2A x iV) = y + y + y “ ax2y ~ yiX ~ y2Xy ’ with “ - t 1 / 4 ) 1 > 
so that the vector field 0 becomes 
/ r \ x = -x 2 + 2 axy + n i + /x 2 y, 

U) o o , 

y = ax — y — y + y2%- 

Elementary calculations and linearization at the fixed points show that, for /2i = 
/i2 = 0, the saddle node remains at (0,0) and a hyperbolic saddle lies at (0,-1). 
Moreover, the y -axis is an invariant line, because x = 0 for any solution with 
initial condition (0,2/o)• The unstable manifold of the saddle (0,-1) is the line 
segment {x = 0 |y G (—oo, 0 )}, the upper part of which coincides with the the 
lower part of the strong stable manifold {x = 0 |y G (— 1 ,+oc)} of the saddle- 
node. A disk containing these two fixed points constitutes a chart, containing the 
codimension 2 degenerate vector field, that can be mapped onto a an open set 
of § 2 : see Figure (8^b). The term —ax 2 y is necessary to make the lower saddle 
hyperbolic (its eigenvalues are —2a and +1). A second hyperbolic saddle lies at 
( 2 a/( 4 a 3 — 1 ), l/( 4 a 3 — 1 ), but this fixed point is irrelevant to the bifurcations of 

interest, and it can be driven out of any compact region by letting a -A ( 1 / 4 ) 1 / 3 ^= f 
a* « 0 . 62996 . For the cases shown in Figure (8^b-j) we set a = 0 . 62996 . 

We now describe the codimension 1 bifurcations and structurally stable vector 
fields that emerge from the codimension 2 bifurcation point for small fi i, /i2- Setting 
x = y = 0 in ([ 5 ]), and noting that y = ( x 2 — [i\)j( 2 ax + /12) from the first equation, 
we may eliminate y from the second equation to obtain the fixed point condition 

(6) G i,iu, 2 ,a( x ) = a 4X 4 + a 3 x 3 + a 2 x 2 + aix + a 0 = 0, 

where 

<24 = 4 <a 3 — 1 , <23 = 2a/4a/i2 — 1 ), cl 2 = 2 /jli + ba/j^ — fi 2 ^ 

( 7 ) <21 = 2 afi\ + fi2 and <20 = /21/22 — 

For fi\ = 112 = 0 Eqn. (| 7 |) becomes (( 4 a 3 — l)x — 2 a)x 3 = 0 , with a triply-degenerate 
root at x = 0 and the irrelevant root at x = 2 a/ ( 4 a 3 — 1). Setting a = a* so that 
the latter root lies at 00, the quartic polynomial becomes a cubic with discriminant 

(8) A = 18aoai<22<23 — 4aoa2 + a 2 <22 — 4a 3 <23 — 27<2o<2 2 . 


To obtain an explicit approximation for the saddle-node bifurcation curve /ii = 
/i(/i 2 ;<^*), we consider this special case. Substituting the expressions 0 into 0 
and setting A = 0 yields a polynomial relating /21 and fi 2 for which (xo) = 

F{n,fi 2 ,a*( x 0 ) = 0 an d one °f relevant roots xq is multiple. Except for fii = fi 2 — 0, 
for which xq = 0 and F" ± a *(Q) = 0, this is a double root, and it corresponds 

either to a saddle-node bifurcation, or to the heteroclinic saddle-saddle connection 
discussed below. Expanding fii in integer powers of fi 2 and using the fact that 
a 3 = —2a + 0 {fi 2 ) to determine the leading terms, we find the following expression 
for the saddle-node bifurcation: 


l4 3 • 2 1 / 3 G 
4 8 


5 • 2 2 / 3 /4 

8 


<9(/4)- 


( 9 ) 


Mi = /i(m2;«*) 
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(c) 




Figure 8. Bifurcations of the gradient vector field <§• (a) The 
bifurcation set in (/ii, /i 2 )-space near the codimension 2 point (0, 0): 
saddle-nodes occur on the curve fi i = fi(/i 2 ', <a) ^ 0; saddle-saddle 
connections exist on the curve Hi = f 2 ( 1 ^ 2 ] = {/^i/4ev 2 |/X 2 €= 

(0, a/g)}. (b) At Hi = fi 2 = 0 the unstable manifold of the saddle 
at (0,-1) lies in the strong stable manifold of the saddle-node at 
(0,0). (c) For /ii > max{/i(/i 2 ; g),/ 2 (/i 2 ; <^)} two unconnected 
hyperbolic saddles coexist with a sink, (d) Along fii = / 2 (/i 2 ; g) < 
a 2 with /i 2 > 0 a codimension 1 saddle-saddle connection exists on 
the invariant line x = (e) For /i(/i 2 ;a) < < $ 2 (^ 2 ] 

and /i 2 > 0 two hyperbolic saddles and a sink exist with the lower 
saddle’s unstable manifold passing left of the upper saddle’s stable 
manifold, (f) On /ii = f 1 (^ 2 ; a) < 0 with /i 2 > 0 a codimension 
1 saddle-node occurs with the lower saddle’s unstable passing to 
its left, (g) For /ii < /i(/i2;o0 only the lower saddle exists, (h) 
On /ii = /i(/i 2 ; a) < 0 with fi 2 < 0 a codimension 1 saddle-node 
occurs with the lower saddle’s unstable entering from its right, (i) 
Along fii = / 2 (/x 2 ;a) with /i 2 < 0 the lower saddle’s unstable 
manifold lies on the invariant line x = — y/JT[ and intersects the 
strong stable manifold of the sink; this is not a bifurcation point. 
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As in Eqn. ©, fi i primarily controls the saddle-node bifurcation, but the second 
parameter y 2 shifts the relative x positions of the upper and lower saddles, allowing 
a codimension 1 heteroclinic connection to form with /1 1 ^ 0. Specifically, along 
the curve 

(10) M 1 = /2(M2;a) = with /x 2 e ( 0 , \/a), 


both saddle points lie on the invariant line x = — yfjli and a connecting orbit from 
the lower to the upper saddle exists (their y coordinates are \ — 1 =p ^/l — [i\joL 

respectively). This bifurcation curve is shown in Fig. |8ja) for a = a*, together 
with the saddle-node curve yi = f 1 ( 112 ; a*) (the latter’s curvature is exaggerated 
for clarity). Note that the discriminant A = 0 for yi = /x|/4cv* 2 since both saddles 
have the same x-coordinate. A similar invariant line x = connects the lower 

saddle to the strong stable manifold of the sink for fi 1 = with /12 < 0, 

but since saddle-sink connections are structurally stable, no bifurcation occurs here 
(Figure gj)). 


In addition to the degenerate codimension 2 vector field at (/ii, H 2 ) = (0, 0) shown 
in panel (b), panels (c-h) show representative vector fields on the codimension 1 
bifurcation curves and structurally-stable vector fields in the three open regions in 
the bifurcation set of panel (a). For the unfolding parameters used here, the saddle- 
node and saddle-saddle connection bifurcation curves meet in a quadratic tangency 
at (/ii,/i 2 ) = (0,0). The geometrical parameters (d, 0) chosen in the construction 
that follows produce bifurcation curves that meet transversely at (0,0), as shown 
in Figure [9j 


4. Geometrical part 


In this section we prove Theorem [3J To do this, it suffices to create for any pri¬ 
mary or secondary edge E = {v\, V 2 } of the metagraph Q a suitable, one-parameter 
family K( A) of convex bodies, where A G [Ai, A 2 ], with a unique value A* G (Ai, A 2 ) 
such that the graph of K( A) is homeomorphic to v\ for any A G [Ai, A*), homeomor- 
phic to V 2 for any A G (A*,A 2 ], and to the graph of the 1-codimension bifurcation 
defined by E = {iq, ^ 2 } at A = A*. In this case we can choose a re-parametrization 
of this family that will satisfy the topological equivalence condition of the theorem. 


We prove the assertion only for secondary edges of Q, because for primary edges 
we may apply a simpler version of the same argument. As noted before, our ar¬ 
gument does not apply to tertiary edges. Secondary edges correspond to non-local 
bifurcations, so it is hard to construct by local truncations a suitable one-parameter 
family of convex bodies that corresponds to any given secondary edge. To ensure 
that local truncations suffice, we rely on Lemma [l] stating that any secondary edge 
belongs to a triangular face of Q of which the two other edges are primary. Since the 
latter correspond to local saddle-node bifurcations, we can use local truncations. 
We will show that any face of Q spanned by two primary edges and one secondary 
edge can be realized by a suitable 2-parameter family K(d, 0) of convex bodies (cf. 
Definition 10). Such a family has (among others) the property that it collapses to 
family described above if we restrict to any of the three edges of Q, so the existence 
of this suitable 2-parameter family proves the Theorem. 
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Figure 9. A codimension 2 bifurcation on a triangular face of 
the metagraph Q. The generic graphs A, B,C can be regarded as 
subgraphs of a topology graph. For example, by adding one sta¬ 
ble point, they are identical to triangulated representations of the 
graphs (f,j,k) in figure [3j The degenerate graphs AB , AC , con¬ 
taining codimension 1 saddle-node bifurcations SN, correspond to 
primary edges fj, fk of Q. The degenerate graph BC , containing 
a codimension 1 saddle-saddle connection H2-H3 corresponds to 
the secondary edge jk of the metagraph. Finally, the degenerate 
graph ABC , containing the codimension 2 bifurcation, corresponds 
to the triangular face fjk of the metagraph Q 


Let B and C be the primary graph representations of two gradient vector fields 
that are connected via any given saddle-saddle bifurcation. Furthermore, let A be 
their common ancestor, that is, B and C can be derived from A by twin vertex 
splittings. By Lemma [l] such a graph exists and from [8] we know that each of 
the three graphs A, B, C can be associated with the gradient vector fields of the 
smooth, convex bodies 104, Kb , Kq , respectively. We denote the degenerate graphs 
belonging to the corresponding transitions by AB , AC , BC and ABC , respectively. 
See Figure [9j 

Definition 10. A 2-parameter family K(d, 0) of convex bodies, where d G [0, dsc] 
and 0 G [0b, 0c] is called suitable if the function (d, 0) i-A K(d, 0) is continuous 
with respect to Hausdorff distance, and there is a value <04 G (0, dsc) an d a function 
0bc : [dA^dsc] £= (0b,0c) such that the following holds: 
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1) 

2 ) 


© 3 ) 

@4) 


5) 

6 ) 
■ 7 ) 
8 ) 


(10 

(10 

(10 

(10 


(10 

(10 


for every 0 G [0b, 0c], K(0, 0) = K A , 

for every 0 G [0b, 0c] and d < <04, the graph of K(d , 0) is homeomorphic 
to A, 

for every d > dA and 0 < 0Bc(d), the graph of K(d, 0) is homeomorphic 
to F>, 

for every d > dA and 0 > 0Bc(d), the graph of K(d, 0) is homeomorphic 
to (7, 

for every 0 < 0bc(^a), the graph of K(dA , 0) is homeomorphic to AF>, 
for every 0 > 0bc(^a), the graph of K(dA , 0) is homeomorphic to AC, 
for every d > <04, the graph of K(d, 0Bc(d)) is homeomorphic to F>(7, 
the graph of iL(d^, 0 bc(^a)) is homeomorphic to ABC. 


If the same properties hold with £de center of mass of Ka as a fixed reference 
point , we say that K(d, 0) is weakly suitable. 


This definition is illustrated in Figure [9| We remark that, in the context of 
Section[3j the line {d = d^|0 G [0b, 0c]} and the curve {0 = 0Bc(d)|d G [d^, dBc]} 
form the bifurcation set associated with the gradient vector field. 

We prove the assertion in two steps: in the first step (Subsection |4.1| ), we con¬ 
struct a weakly suitable family. In the second step (Subsection |4.2| ), we modify the 
construction in such a way that the center of mass of every member of the family 
coincides with that of Ka , showing that the previously constructed family is not 
only weakly suitable but can be made suitable. 


4.1. Neglecting the motion of the center of mass. In the first step of the 
proof we assume that the graph of every convex body is taken with respect to the 
center of mass of Ka , he. we assume that the displacement of the center of mass 
does not influence the topology of the flow. For brevity we set K = Ka , and we 
consider only the case that the equilibrium point of K to be split is stable; if it 
is unstable, a similar argument with an arbitrarily small, conical extension of the 
surface can be applied. 

Let s denote this stable point and the descendant points in the graphs B and (7, 
obtained by splitting s, be s ' B , s B , s ' c , s'f ,, respectively. Appealing to Lemma 5 of 
|5], we may assume that a neighborhood of s in dK belongs to a sphere §. Without 
loss of generality, let the origin o be the center of this sphere, where the radius of § 
is assumed to be one. Furthermore, let c denote the center of mass of iL, and note 
that, because s is a stable point, c is contained in the interior of the segment [o, s]. 

Let r i5 where i = 0,1, 2,..., ra denote the edges of A starting at s, in counter¬ 
clockwise order around s, from outside K. Clearly, for each value of i, the part of 
in § is a great circle arc. These edges are labeled in such a way that the edges of B 
starting at s' B correspond to the Yf s with i = 1, 2,..., k (and those starting at s B 
correspond to the remaining edges), and the edges starting at s' c correspond to the 
T^s with i = 1, 2,..., fc + 1 (and those starting at s'f, correspond to the remaining 
ones). Observe that, measured in counterclockwise order, either the angle from Y i 
to Tfc+i, or the angle from to Y m is less than i r. Without loss of generality, 
we may assume that the angle from Yi to is less than i r. 

First, we truncate the spherical neighborhood of s by a plane P sufficiently close 
to but outside 5, and investigate the equilibrium points of the truncated body with 
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respect to c. In the generic case we have two possibilities for the graph of the 
truncated body Kp. If K D P does not contain a new stable point, then the graph 
of Kp remains homeomorphic to A. Furthermore, if KDP does contain a new stable 
point s", then a new saddle point is created on PndK , and every heteroclinic orbit 
on K intersecting P ends up at s (cf. Figure [To|), whereas those not intersecting it 
remain the same. Finally, note also that K D P contains a stable point if, and only 
if, the orthogonal projection of c onto P is contained in the interior of K D P. (If 
the projection is contained on the boundary of this circle, it is a degenerate case 
corresponding to a saddle-node bifurcation.) We will find a 2-parameter family 
of planes such that, if the intersection circle contains the projection of o on any 
member, then the edges meeting the circle are either Ti,..., Tfc, or Ti,..., Fk+i: 
see Figure [T0| 



r* 


Figure 10. Truncation by the plane P(dpc , 0, </>) to create a body 
with graph B containing a new saddle on P D dK and a sink s ". 

An arbitrary plane in 3-space, and thus, in particular, the truncating plane, can 
be defined with three parameters. For this purpose we use the following coordinates: 

(i) d : the depth of the cut (i.e. the height of the truncated spherical cap), 
measured from the point of the sphere where the tangent plane is parallel 
to the cutting plane; 

(ii) 0 : the arc distance of the center cp of the intersection circle (the one created 
by the cutting plane on the sphere), from s, measured on S; 

(iii) (j): the angle of the great circle arc between cp and s, and from some fixed 
great circle arc starting at s. 

Up to a linear transformation, these parameters correspond to the polar coor¬ 
dinates of the vector pointing from the origin o to its orthogonal projection onto 
the truncating plane P, where the North Pole of § is s. Henceforth we denote the 
plane by P(d, 0, (j)). 

Observe that, measured in counterclockwise order, we have Z(Ti, T^) < Z(T i, Tk+i 
Z(r 0 ,r fc+ i) < Z(r 0 ,r fc+2 ). Choose some angle 0 < a < tt satisfying 

(ii) z(ri,r*+i) < a < z(r 0 ,r fe+ i). 

Furthermore, for any sufficiently small, fixed value 6 > 0, there is a value dpc = 
d B c(0, a) independent of (j) such that for any plane P with parameters P(dpc , 0), 
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a is the angle between the two great circle arcs on the sphere, starting at s and 
touching the intersection circle. Hence, by (11), there are some pB < Pbc < Pc , 
with pB and pc depending only on <a, and Pbc = pBc(d) depending on a and d, 
such that 


• for any p G [pB,pBc(d)) the plane P(cIbc,0,P) intersects if and 
only if i m 1 , 2, ..., k; 

(12) • the plane P(d#c, 0, pBc(d )) intersects if and only if i = 1, 2 ,..., fc, 

and it is tangent to 

• for any p G (Pbc (d), Pc\ the plane P(cIbc,0,P) intersects if and 
only if i s= 1 , 2 , ..., k + 1 (cf. Figure [To]) . 


Now, consider the one-parameter family P(dsc , 0, 0), with # fixed and depending 
only on p G [0s, 0c]- If, for any value of p in this interval, the projection of c lies 
on P fl P(dsc,^,0), then, depending on the value of 0, the graph of the body 
truncated by the plane is homeomorphic to either B or (7, or in the degenerate case 
to BC (cf. Figure [To|). Since we intend to use local truncations only, we would 
like to guarantee this property for any sufficiently small value of 6 > 0. Before 
proceeding further, we recall two lemmas from [B . 


Lemma 2. Let r > \s — c\ and 5 > 0 be arbitrary. Then there is a convex body 
K' C K satisfying the following: 

(i) The graph of K' is homeomorphic to A. 

(ii) Denoting the critical point of K' corresponding to s by s', s' has a spherical 
cap neighborhood in dK', of radius arbitrarily close to r. 

(iii) Denoting the integral curve of K' corresponding toTi by T[ for every i, and 
by ti and t\ the unit tangent vectors of Vi and T[ at s', respectively, we have 
that 11 \ —ti\ <5. 


We note that the same statement is proven in [8] for the case that s is an unstable 
point, and the radius of its spherical neighborhood is arbitrarily close to any given 
value 0 < r < \s — c\. 

Lemma 3. Let C be the unit circle in the plane M 2 with the origin o as its center, 
and let c = (0, r), where r > 0. Let q\ = (fi i, z/i) and = (/i 2 , ^ 2 ) be two points of 
C such that Vi > 0. 

(i) If [qi , ^2] is perpendicular to [s,qi\, then lim — = 

H\ —A) 1 r 

(ii) If the angle of [qi, <72] an d [c, qi] is if-C/ii for some constant C' independent 

of /ii, then lim ^ + 2(7'. 

^1—>-0 ^ 1 1 r 


Now, consider a plane P = P(dsc,^,0) with an arbitrary value of p (cf. Fig¬ 
ure 10), and let the closest and the farthest points of the circle C = PndK from the 
segment [o, s] be denoted by q\ 2 and q\, respectively. For convenience, we imagine, 
for the moment, the plane containing q\ — (/ii,z+), q\ 2 = (/a 2 ,^ 2 ) and c = (0,r) as 
M 2 in Lemma [ 3 ] Figure 11 


For any sufficiently small 0 > 0, we require that the orthogonal projection of c 
on P lie in the interior of the segment [gi, # 2 ]- Since Z(qi,q 2 , c) < this property 
holds if and only if A(q 2 ,qi,c) < | for any sufficiently small 0 > 0. Recall that 
dsc is defined by the fact that the angle of the two tangent lines of (7, passing 
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Figure 11 . An illustration for Lemma [3j The sinks s and 


(filled circles) and the saddle (circle with cross) of Figure 10 
identified. 


are 


through p , is equal to a. Let (7 s , q{ and denote the central projections of (7, 
qi and q ^, respectively, onto the tangent plane of K at s. Then, as 0 -> 0, the 
limit of the angle of the two tangent lines of (7 S , passing through s, is equal to a. 
Thus, an elementary computation yields that, as 0 -A 0, the limit of the ratio of the 
x-coordinate of to that of qf is equal to | , implying that the same holds 
for lim 


We conclude that our requirement that the orthogonal projection of c on P lies 
inside [<71,^2] for any sufficiently small 0 > 0 is satisfied if < i_sin(a/2) ? not 

if d^y > i^sin(a/ 2 ) • To g uaran t ee the former, we apply Lemma 3 and choose d > 0 
sufficiently small, i.e., such that for the truncated body K' and heteroclinic orbits 
T^, the inequalities (11) remain true with the same value of a. Let c' be the center 
of mass of K' and o' be the center of the spherical neighborhood of s. Furthermore, 
let r' = |——— 7 j. Note that for a suitable choice of r, we have > ^-+ sm ( Q /^) 

|s —o | 5 1— r' l-sin(a/2) 


According to the previous paragraph, with a little abuse of notation, we assume 
that for the original body K , for any sufficiently small 0 > 0 , the orthogonal pro¬ 
jection of c on P = P{d BC , 6 , 4>) lies in the interior of P D K. Let = d^(d, a) 
denote the value of d, independent of 0, at which the projection of c lies on the 
boundary of P fi K. 


We have shown that, for 0 G [0b, 4>c] an d 0 is sufficiently small, the intersection 
circle P(d B c , 0) FI K contains in its interior a new stable point with respect to 

c. Thus, the graph of any such truncated body K is homeomorphic to either B or 
(7, or to BC. Then we fix a sufficiently small value of d, and take the 2-parameter 
family of convex bodies iF(d, 0), where d G [0, d^c], an d 0 G [0b, 0c], defined 
as the truncation of K by the plane P(d, d,0): see Figure 12). Finally, for any 
value of 0, iF(O,0) = K , which shows that (fIo}l) in Definition 10 is satisfied. The 
remaining properties in Definition [lO] of a weakly suitable family follow from (3). 
This completes the first step of the proof. 


We note that the bifurcation diagram of Figure [9] in the geometric parameters 
d, 0 used for the construction of the truncating plane and that of Figure [SJa) in 
the unfolding parameters /ii,/i 2 of Section [3] are topologically but not differentiably 
equivalent. As noted in Section [3J the bifurcation curves meet in a tangency in 
Figure [8j[a) ; however, they meet at a nonzero angle in Figure [9j 
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Tangent planes, corresponding to d=0 


Figure 12. The 2-parameter family of truncations used in the 
construction. Variation of the angle of rotation of the truncating 
plane results either in a saddle-saddle bifurcation, or in no bifur¬ 
cation. Variation of the depth d of the truncation results in a 
saddle-node bifurcation; the graphs belonging to the two extremal 
values of <p are identified with capital letters. 


4.2. Annihilating the motion of center of mass by an auxilary truncation. 

In this subsection we modify the family K{d,(j)) in such a way that the center of 
mass of every member in the modified family remains at c. To do this we need 
some additional assumptions on K. 

Let L be the line passing through s and c, and let w denote the point of L D dK 
different from s. We show that K = Ka can be chosen in such a way that q is 
not an equilibrium point, and that it does not belong to any edge of A. First we 
modify the convex body Kq in class (1,1) in [28j to satisfy this property. Since the 
graph of Kq does not contain edges, we need only show that no line through the 
center of mass passes through more than one equilibrium point. 

Since Kq has D4 rotational symmetry, in a suitable coordinate system, its two 
equilibrium points and center of mass c he on the z- axis, and Kq is symmetric with 
respect to the (x, ^-coordinate plane. Thus, all the tangent planes of Kq, parallel 
to the x-axis (i.e. satisfying the property that one of their translates contains the 
x-axis), touch Kq at points in the (y,z)~ plane. Clearly, cutting off a sufficiently 
small part of Kq near the positive half of the x-axis does not change the number of 
equilibria nor the primary equilibrium class {1,1} of the body. The center of mass 
d of the modified body K' 0 is in the open half space {x < 0}. Hence, if the tangent 
plane of Kq at some point p is perpendicular to the segment [c' , p] , then the outer 
normal vectors of this plane have positive x-coordinates, implying that p is in the 
open half space {x > 0}. To show that any graph A can be associated to a convex 
body Ka satisfying this property, we observe that, by [8| Theorem 1], Ka can be 
obtained from Kq by a finite sequence of local deformations. 






22 


G. DOMOKOS, Z. P.J. HOLMES AND Z. LANGI 


In [8] we also showed that a neighborhood of any point of a non-isolated hetero- 
clinic orbit, or a sink, or a source can be truncated by a sphere without changing 
the class of the graph of the body. Furthermore, by [8j Lemma 1], we obtain that, 
applying a sufficiently small truncation at re, the line connecting the modified sta¬ 
ble point and the modified center of mass intersects this spherical surface. Thus, 
we may also assume that a neighborhood of re is a sphere S'. Nevertheless, note 
that the center of S' is not necessarily on the line L. Let denote the center 
of mass of the truncated spherical cap G(d, (j>) near s. To obtain a modified body 
K f (d, </>), we truncate Ka near q by a second plane P'(d', 0', 0'), such that the cen¬ 
ter of mass of the union of P, and the second truncated (open) spherical cap G', is 
c (cf. Figure [l3|. Clearly, in this case the center of mass of the doubly truncated 
body K'(d, (j)) is identical to the center of mass c of K. 



Figure 13. The second truncation near the critical point w op¬ 
posite to s; the circular arcs he on the spherical caps G and G'. 

Let L(d, (/)) denote the line connecting c and c(d, (j)). First, let d' be fixed. Then, 
changing 6 ' and 0', the locus of the centers of mass of G' is a part of a sphere S^,, 
concentric to S', and the radius of this sphere depends on d' and S' only. Thus, if 
6 > 0 is sufficiently small, for every line L(d, </>) and every (small) value of d! there 
is a unique position of G' such that its center of mass lies on L(d, 0). Let us call 
this cap G' = G'(d'). Note that the center of mass of the union of G(d, 0) and 
G'(d') is c if and only if, the torques about c exerted by the two caps are equal. 
Here, the distance of the center of mass of G'(d') from c is approximately \q ■*— c|; 
that is a fixed value. Thus, by continuity, for every pair of values d, </>, there is 
at least one value of d! such that the center of mass of G(d, (j)) U G'(d') is c. Let 
G'(d, 0) be the spherical cap G'(d'), where d! is the smallest value for which this 
property holds. Then, clearly, G'(d, 0) depends continuously on d and 0, and the 
2-parameter family K \ ( G(d , (j)) U G'(d, 0)) has the required properties. 

5. Summary 

In this paper we showed that the secondary classification of smooth convex solids, 
based on the Morse-Smale complexes of their gradient vector fields, is not only 
complete in the sense that all combinatorially possible Morse-Smale complexes can 
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be realized on smooth, convex bodies, but it is also complete in the more gen¬ 
eral, ‘dynamical’ sense that all generic transitions between Morse-Smale complexes 
represented by non-isomorphic abstract graphs can be realized on one-parameter 
families of convex bodies. Among trajectories of physical convex shape evolution 
processes we find examples of such transitions, so our result implies that from a 
purely geometrical viewpoint, there is no restriction on these trajectories. 

Theorem [3] admits only one-parameter families exhibiting one single bifurcation. 
However, if we only admit saddle-node bifurcations then, based on our argument in 
Section[4]we can formulate a more general claim. A codimension one, generic saddle- 
node is either a creation or an annihilation , depending on whether the number of 
generic critical points increases or decreases by two. As stated before, at saddle- 
saddle bifurcations the number of generic critical points does not change. 

To formulate the claim we introduce 

Definition 11. A generic, one-parameter family v(A) of gradient vector fields on 
the 2-sphere is called strictly monotone if it contains either only creations or only 
annihilations and it does not contain any saddle-saddle bifurcations. 

Using this concept, we can state the following corollary to Theorem |3j 

Corollary 1 . For any generic, strictly monotone, one-parameter family v(A) of 
gradient vector fields on the 2-sphere there exists a one-parameter family K( A) of 
(not necessarily smooth) convex bodies such that Vr^(A) is topologically equivalent 
to v(A) for every value of A. 

To extend this statement further, we make 

Conjecture 1 . Every equivalence class on the family of convex bodies, defined by 
the tertiary classification system, is connected. That is, for any two convex bodies 
K\ and K<± with the same topology graph A there is a one-parameter family K( A) 
of convex bodies, where A E [0,1], such that K( 0) = K \, K(l) = K 2 , and the graph 
of K{ A) is A for every value of A. 

If Conjecture [l] is true, Corollary [l] can be extended to include not only strictly 
monotone, but also generic families. Although our techniques do not admit the 
investigation of tertiary edges of Q, we also formulate 

Conjecture 2. All tertiary edges of Q are physical. 

Regarding geophysical applications, we remark that in primary class {2,2} one 
of the secondary classes (that of ellipsoids) appears to be dominant and the other 
appears to be entirely missing among natural pebble shapes. Our results show that 
one could continuously transform members of one class into members of the other 
class. Apparently, this process exists in natural abrasion only in one direction. 
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